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ELECTRICAL ANALOGIES AND ELECTRONIC 
COMPUTERS: SURGE AND WATER 
HAMMER PROBLEMS 


the course extensive investigations the fields speed and pressure 
regulation the Massachusetts Institute Technology (M.I.T.), Cam- 
bridge, use has been made certain electrical-fluid analogies. Also, aid 
obtaining direct solutions while retaining the basically nonlinear features 
the components the systems being studied, electronic computer has been 
found useful. 

This paper describes the foundation for these techniques, the types 
analogies, methods developed, and equipment employed, together with few 
representative results and conclusions drawn from these studies. The par- 
ticular cases pressure transients uniform pipe and surges simple tank 
have been selected for discussion. However, such restrictions are inherent 
the methods developed. 

Emphasis has been placed the utilization analog techniques extend, 
refine, and clarify analytical procedures secure more thorough understand- 
ing the basic hydraulic phenomena and furnish results that are generally 
useful and accessible. Brief mention made some improved analytical 
procedures developed connection with these studies. 


INTRODUCTION 


The interconnection hydraulic and steam power plants within modern 
electric power networks gives rise numerous complex problems concerning 
both the influence load fluctuations and frequency control equipment the 
stable operation the generating units, and, conversely, the effects the 
transient behavior the hydraulic and mechanical components the units 


comments are invited for publication; the last discussion should submitted 
February 1953. 

Asst. Hydr. Eng., Dept. Civ. and San. Eng., Massachusetts Technology, Came 
bridge, Mass. 


| 


ELECTRONIC COMPUTERS 


the performance the electrical network. representative hydroelectric 
system shown schematically Fig. which the gate opening, 
the head, the power, the speed, and the discharge. subscripts, 
denotes ‘‘conduit,” denotes line,” denotes and 
denotes the course this paper, use will made dimen- 


(a) PICTORIAL DIAGRAM Load 


Surge Tank 


Turbine 


(b) SCHEMATIC DIAGRAM 
Draft Tube 


Surge Tank 


Fie. INSTALLATION 


sionless incremental variables—the variables familiar electrical 
engineers. These will generally denoted the lower case letters cor- 


which the numerators each case represent changes from refer- 


ence condition and the denominators correspond this reference index con- 
dition, denoted the subscript zero. 

Most earlier investigations the regulation field were not extended much 
beyond the point which immediate and practical questions were answered. 
Thus, the classical researches into water hammer and surge phenomena that 
formed the subject many pioneering remained adequate 


Sciences St. Petersbourg, Vol. 8th Serie, 1898. 

Mecanique, Janvier Mars, 1904, 

del colpo Allievi, Reale Academia dei Anno CCCIX, 1912. 

Surge Tank Water Power Johnson, Transactions, ASME, Vol. 30, 1908, 

ASCE, Vol. 83, 1920, 707. 

Théorie des Wasserschlosses bei selbstatig geregelten Thoma, Olde 
burg, Munchen, Germany, 1910. 
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long the various generating stations were loosely tied and operated 
relative isolation. However, with the growth interconnection the level 
which nearly all the plants particular system and even entire systems are 
coupled together, these earlier methods analysis and special solutions have 
become inadequate. 

Straightforward extension and generalization the pioneer studies are 
hampered the presence inescapable nonlinearities the various physical 
components the problem. These nonlinearities occur the hydraulic 
elements hydroelectric plant, for example, the friction loss the conduit 
and penstock, the flow through the wicket gates and tank orifices, and the 
turbine power and efficiency characteristics. Basically, all physical systems 
are nonlinear and all analyses arising through linearization procedures are but 
approximate, agreeing with fact only for those cases which the deviations 
resulting from nonlinearity are insignificant; the present context such 
usually the case for equilibrium stability determinations and for small changes 
operating 

However, there are many instances which linear treatments are not 
sufficiently accurate, including, for the hydraulic system, surges and pressure 
swings for large changes demand flow. Problems such these require 
graphical numerical solutions, which techniques are used similar those 
employed for water hammer surge tank 
However, these basic methods must annexed expressions that reflect the 
actual plant and load characteristics well. The tediousness and consumption 
time inherent these methods emphasize the need for more effective tech- 
niques for solution, without disregarding essential nonlinearities. 

first step toward developing new techniques and refined solutions 
this field, both for the purpose extending basic knowledge and answer im- 
mediate questions arising practice, the M.I.T. Hydrodynamics Laboratory 
began, 1947, investigation into the basic dynamic phenomena involved 
the transient performance power system prime movers, which necessarily 


Gaden, Editions Concorde, Lausanne, Switzerland, 1945. 

stabilité d’un réglage automatique vitesse par des diagrammes vectoriels,” 
Gaden, Informations Techniques Charmilles No. Geneve, Switzerland, 1946. 

certaines caractéristiques intervenant dans condition stabilité. propos 
réglage automatique vitesse des turbines Gaden, Editions Concorde, Lau- 
sanne, Switzerland, 1949. 

Houille Blanche, Novembre, 1945. 

des phénoménes coup pélier sur réglage vitesse des turbines 
Cuenod, ibid., Mars-Avril, 1949. 

regolazione delle turbine Evangelisti, Zanichelli, Bologna, 1947. 

stabilita regolazione nelle installazioni Evangelisti, Elet- 
trica, 

Scimemi, ibid., 1947. 

graphique générale calcul des propagations d’ondes Bergeron, Memoires 
Societé des Ingenieurs Civils France, 1937. 

Hammer Pipes, Including Those Supplied Centrifugal Pumps,” Angus, Pro- 
ceedings, Inst. Mech. Eng., Vol. 136, 1937, 245. 

Schnyder, Schweizerische Bauzeitung, Vol. 94, Nos. 
and 23, 1929. 

chambres d’équilibre,’’ Calame and Gaden, Gautier-Villars, Paris, France, and 
Concorde, Lausanne, Switzerland, 1926. 

Surge Chamber Hydroelectric Cole, Selected Engineering Papers, 
No. 55, Inst. London, England, 1927. 

Berechnung von Braun, Schweizerische Bauzeitung, Band 86, 1925. 
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required research into the field hydraulic transients. Certain these studies 
are outlined here well related paper the author.” 


ANALOGS AND COMPUTERS 


One fruitful approach toward understanding these problems has been found 
the use electrical analogs and computers, means which solutions may 
obtained rapidly and immediately useful form. These devices fall natu- 


rally into two analogs and computers. 
analogs the prototype system, effect, duplicated 


model system whose equivalent components behave according laws analog- 
ous those governing the prototype. 

computers the basic algebraic and differential equations 
are solved either (a) Digital computation (numerical calculation dis- 
crete steps) analog computation (involving calculation continuous 
variation analogous variables). 

Comparisons.—Generally, then, the establishment formal analogy be- 
tween the basic equations describing the behavior two different physical 
systems permits the previously determined solutions for new cases. 
Indeed, this one the major benefits derived from mathematical analysis 
itself. 

problems become more complex, with more variables and system 
parameters, the utility analog techniques becomes more marked. par- 

ticular, most mechanical (or dynamical) problems can represented suit- 
able electrical analogs. Since electrical terminology and concepts have be- 
come highly developed the fields unsteady and periodic motion, these 
analogies will usually prove fruitful. For detailed introduction the basic 
structure analog techniques the reader referred the 

With respect computers, some form high-speed computer useful 
the system under study involves large number variables many different 
solutions for varying parameters and conditions. Even linear analyses became 
awkward for these cases. Computers are also helpful the system possesses 
one more significant nonlinear features that render analytic solutions 
impossible. 

Many problems the field speed and pressure regulation possess both 
these attributes. this paper the latter aspect has, been emphasized, 
especially since number novel results have been obtained. 

and Results from M.I.T. Studies Unsteady Flow,” Paynter, Journal, Boston 
Soc. Civ. Engrs., XIX, No. April, 1952. 


York, Y., 1950 
“Calculating Instruments and Hartree, Univ. Illinois Press, Urbana, 
Mathematical Murray, Kings Crown Press, New York, Y., 1947. 
York, Y., 1942. 
Mathematics for Engineers and Physicists,” Pipes, McGraw-Hill Book Co., Inc., 
New York, Y., 1946, Chapter VIII. 
Inc., New York, Y., 1940, Chapter VI. 
Olson, Van Nostrand, Inc., New York, Y., 1943. 
Glenn Murphy, Ronald Press, New York, Y., 1950. 


Speed Computing Devices,’’ Eng. Research Associates, McGraw-Hill Book Co., Inc., 
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CONCEPTS 


All continuous devices transport energy over distances possess many 
features common. Such mechanisms may mechanical shafts, electrical 
transmission lines, pressure pipe lines, open channels, all which can 
conceived series inertial elements linked together elastic flexible 
couplings capable storing potential energy. other words, these systems 
may visualized terms simple dynamical model consisting mass cars 


and coupling springs shown Fig. energy terms associated 
with these elements are 

and 


which kinetic energy and potential energy when the cars are 
moving with velocity and the springs are compressed distance 

the leading car suddenly stopped, the springs must each turn absorb 
the kinetic energy the adjacent car, and convert into potential energy. 
Furthermore, definite time delay occurs each unit, order decelerate 
the mass and build the compression. These two features may obtained 
from the conservation energy and momentum principles follows: 


Energy— 


(2b) 
Momentum— 
From which 
and 
(6) 


—> V. 
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which the propagation velocity. the case elastic pressure waves, 
the corresponding values for the surge impedance and propagation velocity 
become 


which the acceleration due gravity, and for unbounded fluid 
Ey ‘ 
(8) 


which the density, with relating the change head the change 
velocity. These concepts prove generally useful for all types transmission 
problems. 


ELEcTRICAL ANALOGY 


demonstrated Table that uniform, frictionless pipe line formally 
analogous uniform, dissipationless electrical transmission line. Thus the 


TABLE ANALOGY 


Line (a) Hydraulic System: (b) Electrical System: 
Uniform frictionless pipe line Uniform lossless transmission line 


Wave equations: Wave equations: 
Refiections: Reflections: 
Open end: Pressure node: Grounded end: Voltage node: 
Closed end: Velocity: Open end: Current node: 
ANALOGY 
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water hammer waves and surges the hydraulic engineer became the travelling 
waves, electrical surges, and switching transients the electrical engineer. 
The fruits this analogy are (1) the ability make use the many useful 
tools and concepts developed the electrical engineer since 1900 and (2) the 
realization that, when viewed this light, hydraulic and electrical engineers 
have similar problems and common language. 

voltage 

long dissipation phenomena are negligible both systems these 
analogies are strictly valid. Nevertheless, the practical case resistance and 
losses must taken into account. For the complex forms pipe networks, 
even steady flow, computers are useful. With steady flow one successful 
method making friction proportional exponential power the flow 
embodied the pipe flow For unsteady flow problems, using elec- 
tronic computers, possible account for this frictional effect satisfactorily 
through special components described subsequently. 

However, even the loss-free form, many concepts electrical engineering 
have direct application problems unsteady flow. For example, the 
ordinary alternating-current vector diagrams furnish valuable clues the 
behavior pipe lines subjected alternations flow and pressure. Moreover, 
the use surge impedance and the other generalized circuit constants can 
extremely profitable. 


RESONANCE PHENOMENA 


excellent example the use electrical concepts for unsteady flow 
problems can shown connection with the resonance phenomena asso- 
ciated with the rhythmic motion gate valve the end pipe line. 
such cases, the resultant pressure fluctuations along the pipe may often sub- 
stantially excess the maximum the valve for the full gate stroke, for 
which the governor gate mechanism timing usually specified and adjusted. 
will also interest compare the results obtained these techniques 
with those obtained following the classical methods. 

The solutions the electrical wave equations for long alternating-current 
transmission lines, derived any standard work power 


may put the form: 


and 
0 


which the distance along the transmission line from the sending end; 


the voltage vector point the voltage vector the sending end; 


Electrical Analogy for Pipe Malcolm Mellroy, Proceedings-Separate 
No. ASCE, 1952. 
Electric Power Woodruff, John Wiley Sons, Inc., New York, 

System Rissik, Pitman, London, England, 1940. 

“Standard Handbook for Electrical edited Knowlton, McGraw-Hill Book Co., 
Inc., New York, Y., 1941, Sect. 13. 
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the current vector point the current vector sending end; 

the vector propagation constant; and the vector surge impedance. 
the line short-circuited (analogous open reservoir) then 

and the sending end voltage and current are related the expression, 


(11) 


Furthermore, the line resistance negligible (analogous the frictionless 
pipe), the propagation constant given 


which the angular frequency transmission; the propagation ve- 

the analogy between flow velocity current, and head voltage, 
one may rewrite this expression directly hydraulic terms, obtain cor- 
responding values for head and velocity the gate, thus: 


and 


which equals natural period pipe; the period gate oscil- 


lation; equals the Allievi parameter; equals the relative 
head change; and numerically equal equals the relative ve- 


locity (or flow) change. 

Eqs. and demonstrate that the head variation will always lag behind 
the flow variation angle 90° and that the amplitude ratio head rela- 
tive flow given the expression: 


For small variations, however, the normalized flow may also found from 
the expression: 


which relative gate opening, or, approximately 
Combining Eqs. and 17, 
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and 


graphical representation these expressions shown Figs. and 
These results were first obtained using intervalfequations derived 
Allievi. 


Values 


Values of K= 2oHy 


RESONANCE 


For example, Fig. assume =.5 and 2.5 locate the point 


then extend line through point point The values 
1.92 and 164° are read directly. The significance these curves 


arises from consideration the variations velocity (or flow) and input 
torque (or power) caused oscillations gate opening. The velocity and 


Phase 
Period T, 
Time 
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input torque the normalized (or per-unit) dimensionless form may ex- 
pressed the equations: 


Flow— 
and power— 
Values 


" 
90° 
DIAGRAM 
For small increments Eqs. and become 
and 


These sinusoidally varying increments can conceived rotating vectors (as 
illustrated Fig. and the resulting amplitude ratios, making use the 
previously determined expressions, become 


and 


with the lagging angle between the gate vector and the head vector 
before, increasing from 90° (for 180° (for 1). 
Thus, both the velocity and the power always lag behind the gate 
demonstrated Fig. which should noted that the tips the v-vector 
and p-vector, well the multiples all execute circular loci the fre- 
quency increases. 
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(9) 
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may concluded from these analyses that the input power the prime 
mover inherently lags behind the gate opening, matter how slowly the gate 
moved, with increasing phase shift for higher frequencies; this circum- 
stance critical determining the transient performance hydroelectric 
units and clearly indicates one the fundamental limitations such units for 
participation system load and frequency control programs. 

the Appendix, numerical example given showing the use Figs. 
and and comparing the results with those obtained conventional methods. 


APPLICATION RESULTS 


From the foregoing evident that oscillations the gate valve produce 
standing waves the pipe systems pipes with the familiar loops and nodes 
all vibratory continuous systems. the gate end there will usually 
pressure loop, and the upstream end there node the free surface. 

From purely intuitive considerations clear that only resonances the 
fundamental and odd can occur under such conditions; but for all 
these cases, the full pressure variation will occur each loop, making the 
resonance the harmonics, general, more dangerous than that the funda- 
mental. 

Moreover, can shown that possible, under certain conditions, 
produce pressure swings along the pipe (through oscillations the gate that 
are consistent with the governor timing and stroke) nearly large the design 
values established the conventional methods. These extreme fluctuations 
can exist both for the fundamental and the lower harmonics. For example, the 
pipe line illustrated Fig. 6(a) would serious danger bursting col- 
lapse near point under the conditions shown. 

Another dangerous case simultaneous resonance can arise between the 
gate, penstock, and conduit shown Fig. this case, the surge 
tank will not serve trap the elastic waves since node exists its base, and, 
therefore, renders the tank inoperative. Failure resulting from excessive 
positive negative pressures may occur the loops. 


Considerable irsight analysis and behavior simple surge tanks 
can gained through consideration the energy principles involved the 
elementary problem sudden, full-load rejection with cylindrical tank. 
The basic dimensions and physical constants such tank are indicated 
Fig. 

For the limiting case frictionless conduit, the conservation energy 
may expressed the form: 


or, specifically, 


with and being the mass water the conduit and tank, respectively. 
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With the data Fig. there results 


L.) y2 (w A: 24, Le) V% (28) 


which and are constants. 


c Gradient 


Dynami 


Tank 


This expression therefore the equation ellipse which can always 
transformed into circle the process normalization that is, making 
horizontal and vertical scales homogeneous. Thus, dividing through the 
term 


a V% a V% 


which 


with defined the surge impedance the tank. This situation 
t 


sketched Fig. 8(a), which seen that maximum positive surge 
occurs point for which 


From Fig. 8(a) also clear that the tank will oscillate indefinitely result 
continuous interchange energy without dissipation. 


L, 
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the practical case, however, conduit friction will modify the basic energy 
expression into the form: 


(work done against friction)............ (35) 


Thus the total change water level the tank will increased the 
presence friction, although the surges will dampen out successively result 


Relative Tank Level Level 


Conduit 
Friction Loss 


WITHOUT FRICTION (6) WITH FRICTION 


The effect friction, for any given transient, can measured the ratio 
friction loss free surge known. Thus the parameter may defined 


(36) 


which the rated conduit friction loss. For particular value the 
surge oscillations would modified shown. should noted that the 
initial friction loss relative given R/2. 

has been found useful and convenient multiply both the horizontal 
(conduit flow and vertical (tank level scales the constant factor 
obtain the normalized plots shown Fig. 

Further analysis will show that terms the notation used the dia- 
gram— 

becomes— 


Eq. resolves the problem computing surges into the evaluation the work 
integral the right-hand side, corresponding the areas and Fig. 
This can done several different procedures and has been successfully 
achieved using the electronic computer. 
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The value the particular formulation the problem implied Eqs. 
and lies the considerable precision that may derived using even ap- 
proximate methods integration. Employing such techniques, tables 
surges that are accurate four significant figures have been prepared for the 
normal range tank design. Similar procedures have been developed for 
restricted orifice and differential tanks. 


Tank Level Tank Level Tank Level 


Maximum Positive Surge Load Rejection 


Initial Conduit 
Steady 


Flow 


Friction 
Loss 


Final Steady Level 
Steady State Friction Drop 
Acceptance 


Final Steady Level 
Maximum Negative Surge Load Acceptance 


(a) SURGE-TIME CURVE (b) SURGE-FLOW SURGE-FRICTION CURVE 


Initial 


Tank TRANSIENTS 


also interest note that the tank constant serves dimension- 
less similitude index similar the Reynolds, Froude, and Mach 
numbers. other words, two simple surge tanks and will dynamically 
similar 


This conclusion agrees with the tank similitude analyses 
and 


PROBLEMS 


The excellent early studies were essentially limited only 
establishing design constants and dimensions for fixed positive and negative 
surge heights. Mr. Johnson did not attempt the determination the various 
relationships between diverse values the many system constants, but rather 
established working design values for particular choice these constants; 
these restrictions appreciably simplified the analysis performance, leading 
directly the type results that Mr. Johnson sought. 


the Surge Tank Problem Model Experiments,” Gibson, Proceed- 
ings, Inst. E., London, England, 1924-1925 

Gibson, ibid., 1932-1933 

actions, ASME, Vol. 43, 1921, 1177. 
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However, many tanks now installed are often required operate under 
conditions quite different from those selected for design, and the behavior 
under these conditions has sometimes been unfavorable and unexpected. The 
problem determining the response tanks under all reasonable disturbances 
that might encountered has been neglected the past, largely because the 
altogether forbidding length time required even the most elementary in- 
vestigations using conventional methods. 

With these needs view, the research efforts reported this paper have 
been devoted, the surge tank field, the following problems concerning pres- 
ent and proposed tank installations: 


The behavior tank systems for operating conditions other than those 
used for design; and 

The true stability margins for surge tanks subjected specified dis- 
turbances appreciable magnitude. 


has long been realized that the actual transient behavior surge tanks 
considerably influenced the action the turbine governor and other plant 
and load characteristics. However, the magnitude this effect depends the 


Tank Level 


Tank Level Tank Level 


Surge 


B,= Constant 


(a) CONSTANT GATE CONSTANT FLOW (c) CONSTANT POWER 


Constant 


AssuMPTIONS 


comparative values the governor timing and the natural period the tank. 
Investigations the writer have established the limitations and range 
validity three simplifying assumptions plant effects, illustrated 
Fig. 10—constant gate opening, constant discharge, and constant input 
power. 


Constant Gate Opening Constant) the Final Steady-State Value.— 
This the situation most representative unit under manual control the 
operator, and produces the least surge, under otherwise comparable conditions, 
result increased positive damping action. 

Constant Discharge (Q2 Constant) the Final Steady-State 
This was the basis the methods Mr. Johnson and most other investigators 
and produces greater surge than assumption 

Constant Input Power Constant) the Runner the Final Steady- 
State assumes that efficiency changes are small and the turbine 
regulates over time that very short compared the tank period. This 
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assumption gives materially greater surges, and can even lead instability, 
explained subsequently Fig. 14. 

This method, when suitably linearized, was the basis the classical tank 
stability analyses’ that have been modified later Calame and 
Gaden, Evangelisti, and the 


addition the conventional problem sudden step-change 
demand conditions, the electronic computer has been applied successfully the 
examination two other types problems—synchronous resonant load 
changes and true static and dynamic stability margins. 


Synchronous Resonant Load Changes, the Pulsing and Oscillatory 
loading program may result abnormally high surges and 
pressure swings the period the disturbance approaches the natural period 
the tank. Previously, this problem had been investigated only for handful 
cases.*° 

True Static and Dynamic Stability Margins for Tanks Operating Both 
Under Steady-Load Conditions and Under Appreciable Changes Demand.— 
The simplified techniques the linear stability theory, for example, cannot 
account for the known augmentation stability through the differential princi- 
ple, giving rather the same value for the critical tank area for all common types 
tank. Moreover, the engineering concept stability implies far more than 
the various classical definitions mathematics and mechanics. surge tank 
stable real sense, must able reach and maintain effective 
steady operating condition sufficiently short time after any reasonable load 
disturbance. Simply having the oscillations die out time just short 
eternity not adequately meeting these requirements, yet quantitative, gen- 
erally applicable knowledge has been lacking these matters. 


With the conclusion first stage the current research program 
possible give few preliminary examples the many results 
obtained. 


TANK 


For computational purposes, usually convenient return the basic 
differential equations the physical form; for the simple tank, following the 
nomenclature Fig. 


Continuity— 
Aimy = Q> Q. Li (41) 
Acceleration— 
dQ. 
(42) 


Assuming uniform conduit and tank, with square law friction, these may 


stabilité des instalations hydrauliques munies des chambres Calame and 
Gaden, Sehweizerische Bauzeitung, Band 90, 1927. 

Cambridge, Mass., 1949, partial fulfilment the requirements for the degree Master 
Science 


Differential Surge Johnson, Transactions, ASCE, Vol. 78, 1915, 760. 
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normalized the form: 


Continuity— 
Acceleration— 


(a) TANK LEVEL VERSUS TIME (b) CONDUIT FLOW VERSUS TIME 


(c) TANK LEVEL VERSUS CONDUIT TANK LEVEL VERSUS FRICTION 


CoMPUTER SOLUTIONS SuRGE TANK TRANSIENTS 


might take the transient from initial steady state 1.7 final 
steady state that might either full partial rejection zero 
for great number practical instances. This transient indicated Figs. 
and 11. 

Thus, any given plant can specified the value the parameter for 
rated conditions, and the particular transients are specified the initial and 
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The work Milne this involved complete investiga- 
tion numerical methods the system Eqs. and Milne also 
described the general case quadratically damped vibration; but his solu- 
tions were only for step change demand flow, assumption 
his tables use the identical normalized variables outlined this paper and 
hoped that this paper may serve another introduction Mr. Milne’s 
valuable work. 


NONLINEAR FEATURES 


The two principal nonlinearities the simple surge tank equations, out- 
lined, lie the demand flow characteristic and the friction loss term u?/2. 

The Demand Flow Characteristic factor may have any the 
three forms indicated Fig. 10, and may vary with time well because the 
finite governor time. 

The Friction Loss Term term has been assumed 
quadratic since analysis somewhat simpler for this case; but general 
exponent different from may have taken into account. Although 
numerical graphical solutions particular instances may desirable, 
has been found possible, means the electronic analog computer, calcu- 
late complete and universal solutions the surge tank equations for all cases 
normally arising practice. 


CoMPUTER SETUP 


The electronic computer used the M.I.T. Hydrodynamics Laboratory for 
research the regulation field the high-speed analog type, and ideally 
suited this type analysis. The basic computing epoch four milliseconds 
and the system variables are represented analogous voltages functionally 
related through appropriate computer components that solve the fundamental 

The general arrangement the computer illustrated Fig. along with 
operational details representative banks components. particular, Fig. 
12(d) shows the arrangement components and the block diagram for simple tank 
transients following sudden changes demand especial interest 
note the C)-bank and its behavior outlined Fig. 12(c); 
this bank that represents the square-law friction. 

Typical computer solutions, photographed from the oscilloscope screen, 
are shown Fig. 11. These may compared the solution diagrammed 
All problems may set such way that calibration and scaling 
can made directly from the final photographic records. This arrangement 

Vibrations; General Theory Together with Solutions Important Special 
Milne, Univ. Oregon Publications, Eugene, Ore., August, 1923. 

Damped Milne, Univ. Oregon Publications, Eugene, Ore., 
March, 1929. 

Derivatives for Damped Milne, Oregon State College Monographs, 
Corvallis, Ore., December, 1935. 


“The Study Oscillatory Circuits Analog Computer Chang, Lathrop, 
and Rideout, Proceedings, National Electronic Conference, 1950. 


and Philbrick Researches Inc., Boston, Mass., 1950. 
Industrial Laboratories, May, 1952. 
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has been found very successful practice and constitutes another strong argu- 
ment for the normalization procedures that make possible. 


ANALYSIS AND RESULTS 


worthwhile present two examples the many charts that have been 
assembled from the computer studies simple surge tanks. Similar curves 
have been obtained for other types tanks. 

The first specimen, shown Fig. 13, family curves, plotted similar 
the method indicated Mr. from which the maximum minimum 
tank level resulting from any load acceptance rejection transient with any 
simple tank may computed, subject the assumptions constant demand 
flow and square-law friction. This plot may compared with the earlier 
Johnson example its use given the Appendix. 


0.2 0.4 0.6 0.8 1.0 2.0 2.2 2.4 
Values 


Tank 


The second curve, Fig. 14, has with surge tank stability for full-load 
acceptance and compares the conditions the classical criterion’ with actual 
computed solutions assuming constant input power, plotting the required 
minimum value the simple, tank constant for any given value conduit 
friction loss relative static head, The curve seen com- 
posed three sections corresponding the three types instability indicated. 
There are, order—the oscillatory region, the drainage region, and the region 
unstable steady state. 

Region (Small Frictional this range stable response 
(below the curve) indicated damped oscillations and unstable behavior 
evidenced oscillations growing magnitude. 

Drainage Region (Medium Frictional this region satis- 
factory behavior results rapidly damped return full-load gradient condi- 
tions. However, the tank insufficient size, the conduit flow would fail 
accelerate fast the tank was being drained, and, the practical case, the 
unit would evidence instability failure deliver the required power the 
bottom the swing. This might result loss synchronism. 


Creager and Justin, John Wiley Sons, Inc., New York, 
Y., 1950, Chapter XXV. 
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Region Unstable Steady State (High Frictional 1/3, 
the static equilibrium unstable and the slightest disturbance would result 
unstable response. Stating the fact another way, the rated friction loss 
greater than one third the static head stable regulation impossible regard- 
less the presence size the tank. 

Comments.—Of course, all these conclusions are based the assumption 
constant power input the prime mover (case However, the effect the 
finite time taken the governor respond the changes head well 
the effects interconnected generating units usually are such ameliorate 
these conditions. Nevertheless, the curves serve the useful purpose design 
limitations and should noted that the conventional allowance for drooping 
efficiency curves applied the Thoma criterion may still inadequate for the 
larger values 

The use these curves shown the Appendix for actual surge tank 
installation which field tests were made check the design 


CoNCLUSION 


hoped that this outline the type methods used the research pro- 
gram will stimulate new interest the desirability obtaining more 
complete and effective solutions problems the regulation field. Prelim- 
inary results have given ample evidence that much valuable information can 
derived through use appropriate analogies, whereby results other fields (for 
instance, electrical engineering) may profitably interpreted hydraulic 
context. 

addition, believed that adaptation electronic computers surge 
and water hammer studies offers new possibilities for obtaining results gen- 
eral and permanent usefulness problems whose solution impractical 
any other means. 
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APPENDIX. TYPICAL EXAMPLES 


Explanatory Notes—The following two examples have been included pri- 
marily demonstrate the use and limitations the design curves presented 
the paper. particular, the second example, treating surges the simple 
surge tank Tallulah Falls, Ga., was chosen from among the many field test 
results principally because the deviations from computed studies (obtained 
electronic computer any other means) represent the typical range un- 
certainty encountered practice. 


Check Computed Values Surges,’’ Eugene Lauchli, Engineering Record, Vol. 71, 1915, 
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not intended that the very brief exposition methods presented 
considered complete itself, but outlined only demonstrate the value 
generalized curves and formulas for rapid approximate computations. 

The material Example drawn from the curves presented Figs. 
and complete explanation and theoretical treatment water hammer 
resonance have been presented Mr. and Mr. The 
confirmation these results graphical solution follows well-known 

However, Example addition the use the surge diagram Fig. 
13, certain approximate formulas for times and alternate surges have been 
employed that were not described the text the paper. These rough formu- 
las have been found give good results when checked against both the writer’s 
theoretical studies and field tests, and are similar certain expressions pro- 
posed the work Mr. explanation the nature and use 
diagrams similar Fig. given the paper Mr. 

Example Water Hammer hypothetical hydroelectric 
plant with the following physical constants analyzed: 


Factor Value 
Rated velocity Vo, feet per 
Wave velocity feet per second..................... 3220 


From these data the elastic constants and obtained. Thus, 


assumed that gate oscillations exist with period 4.8 sec and 
swing 20% range, the problem concerned with determining the resultant 
Thus, the data, concerning gate oscillations, required for using the resonance 
diagram are: The oscillations are assumed sinusoidal about 100% gate open- 

Referring the resonance diagram (Fig. 4), 161° and h/b 1.90. 
The head oscillations are thus found range over (1.90 38%, the 
period for this case being 4.8 sec, with phase lag 161°. 

The foregoing determination can checked graphically reference 
Fig. 15. Thus, for gate swing 20% and period (as before) 
4.8 sec, the head range found 38% with phase lag 165°. 

Although the resonance diagrams presented were derived the basis 
small changes, the agreement with conventional computations for large changes 
satisfactory for most purposes. However, the upward shift the head 
variation curve should noted, since this assymmetry may reach substantial 
proportions particular cases. 
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Example Surge Tank Tallulah design data and general 
computations Table are self-explanatory. The use these preliminary 
data demonstrated for three cases. 

Case 0.425 and 0.—Enter Fig. 13(a) 0.425 (point A), 
read vertically upward intersect (point B). From this point, read 


(6) TIME PLOT 


(a) GRAPHICAL SOLUTION 


RESONANCE 


horizontally point the scale the right, obtaining +0.365. 
read directly from the chart point 


Thus, summary, 


Relative surge ya, feet above the reservoir surface +0.365 
Relative friction drop feet below reservoir surface 0.091 
Relative level change from the initial level 0.456 


These values can converted physical magnitudes most easily through 
the proportional relationship: 


Level change (in feet) relative change (47) 
Friction drop (in feet) friction 


120% 
Head 
115% 
90% 
85% 
4 P 
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Case ity, incre- Level maxi- Free surges® 
(ft per (ft) (ft per (ft) change (ft) 

sec) sec) (sec) 

(1) (2) (3) (4) (5) (6) (7) (8) 
6.13 5.0 23.6 120 6.13 X3.83 =23.5 =0.425 
3.18 1.5 11.9 3.18 X3.83 =12.2 0.425 X3.18/6.13 =0.222 
0.16 (0.2) —6.65 100 —1.44 X3.83 0.425 X0.16/6.13 =0.011 


The tank this example has rectangular section, ft, with section area 
observed friction measurements are given Cols. and 


Surge impedance (Eqs. and 34) 


342 sec. 


Level change (in feet) friction drop (in feet)....... (48) 
0.456 
Thus, level change 5.0 25.0 ft; and time maximum— 


alternate surges for corresponding surge numbers 


108 


vertically upward intersect From this point, read horizontally 
the scale the right, obtaining +0.207. case the initial 
from Fig. other words, 


Relative surge yz, feet above the reservoir surface +0.207 
Relative friction drop feet below reservoir surface 0.022 
Relative level change from the initial level 0.229 


Again, applying the relationships expressed Eqs. and 48, 
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TABLE SuRGE TIME 
Surge 
num- 
Surge Time Surge Time Surge Time 


0.365 2.740 +20.0 108 0.207 4.831 +11.0 —0.100 —5.5 


+0.667 +170 +0.667 +170 

0.293 3.407 —16.1 278 0.182 5.498 9.7 267 +0.069 +3.8 262 
+0.667 +170 +0.667 +170 

0.245 4.074 +13.4 448 0.162 6.165 8.6 437 —0.073 —4.0 432 
+0.666 +170 +0.666 +170 

0.211 4.740 —11.6 618 0.146 6.831 7.8 607 +0.047 +2.6 602 
+0.667 +170 +0.667 +170 


0.185 5.407 +10.1 788 0.134 7.498 7.1 777 


Time first maximum— 


Alternate surges— 


The remaining steps the computation surge time are entered Table 


cases and 


Relative surge feet above the reservoir surface —0.100 
Relative friction drop fc, feet below reservoir surface 0.000 
Relative level change from the initial level —0.100 


Eqs. and yield: 


Time first maximum— 


Alternate Surges (amplitudes Fig. 16)— 


The computations for the remainder case are shown Table 3(c). 
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values surge heights for the three cases given are plotted 
Fig. 16, for comparison with the actual test records plotted the smooth 
curves. The divergences for the load rejection cases are caused primarily 
the slow rate which the demand flow was reduced. the other hand, for 


+20 


A.—Load Rejection From 
6.13 per Sec 


B.—Load Rejection From 
3.18 per Sec 


From 0.16 
1.6 per Sec 


Time, Seconds 


16.—TRANSIENTS TALLULAH TANK 


the load-acceptance case, the actual curves may expected excess 
the constant demand-flow results because the effort the governor pre- 
serve constant power. 
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